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We investigate effects of pseudo-spin population imbalance on Mott phases in ID trapped two- 
component atomic Fermi gases loaded on optical lattices based on the repulsive Hubbard model 
in harmonic traps. By using the density matrix renormalization group method, we numerically 
calculate density profiles of each component and clarify the pseudo-spin magnetism. Consequently, 
we find that all the features from weakly imbalance to fully polarized cases are well described by 
S — 1/2 antiferromagnetic Heisenberg chain under magnetic field. These results indicate that the 
Mott phases offer experimental stages for studying various interacting spin systems. 

PACS numbers: 03.75.Ss, 71.10.Fd, 74.81.-g, 74.25.Jb 



Recently, effects of population imbalance on interact- 
ing fermion systems have been intensively studied in var- 
ious fields as superconductors, atomic Fermi gases, and 
quantum chromodynamic^. The main reason is recent 
drastic developments of experimental techniques in su- 
perconductors and atomic Fermi gases^. In particular, in 
atomic Fermi gases, one can arbitrarily tune the popu- 
lation imbalance, so that not only the so-called Fulde- 
Ferrell and Larkin-Ovchinikov (FFLO) phase^ with a 
spatially-modulated superfluid order parameter but also 
the Chandrasekhar-Clogston limit^ in a large imbalance 
have been explored. 

In cold atomic gases, besides the tunable imbal- 
ance, the optical lattice and the variable interaction are 
like magic arts for condensed matter physicists^. The 
optically-created periodical potential flexibly builds up 
various playgrounds. The interaction tuning associated 
with the Feshbach resonance provides a chance to sys- 
tematically study strongly-correlated behaviors^. In this 
paper, we therefore study the population imbalance ef- 
fect on the strongly-correlated lattice stage, which is now 
one of the most intensive but controversial issues in solid 
state matters^. 

The atomic gas experiments usually employ the har- 
monic trap produced by magnetic fleld and/or optical 
method to avoid the escape of atoms. The harmonic 
trap brings about spatial inhomogeneities, which compli- 
cate the observation of the quantum phase transition^. 
Moreover, the fact that the most convenient probe is 
atomic density proflle have limited the exploration of 
novel phases^. For example, the sign reversal in the 
FFLO superfluid order-parameter can not be directly rec- 
ognized by the density proflle. Thus, the experimental 
conflrmation of FFLO still remains controversial in the 
trapped systen^^*^. 



On the other hand, the Mott insulator core accompa- 
nied by metallic wings predicted in the trapped optical 
lattice in the presence of the repulsive interaction^ can 
be easily conflrmed by the current probe like the density 
proflle. These inhomogeneous phases have been proposed 
by Quantum Monte Carlo studies^ as well as the exact 
diagonalization method^^. So far, theoretical studies of 
the Mott core phase have been restricted to a particu- 
lar case, "balanced population" . In this paper, we focus 
on the Mott phase in the presence of population imbal- 
ance. Using the density matrix renormalization group 
(DMRG) method^^'^^, we investigate pseudo-spin struc- 
tures by calculating density proflles of each component in 
the Mott phase. Since the Mott core and its pseudo-spin 
structures are directly observable, their exploration will 
be a suitable next challenge in cold atom physics. 

Inside the Mott core, the on-site atomic density shows 
the unit-fllling and the density compressibility vanishes^. 
As a result, the pseudo-spin degree of freedom solely sur- 
vives, so that the core region is well described by 5* = 1/2 
Heisenberg (local pseudo-spin interacting) model for the 
two-component atomic Fermi gas. Moreover, we expect 
that a population imbalance has a role of the magnetic 
fleld in the Heisenberg model given by 



(id) ^ 



(1) 



where the flctitious magnetic fleld i^ext is varied by 
the magnitude of the population imbalance in the orig- 
inal system. In this paper, we explicitly conflrm that 
the spin structure in the Mott core region is really de- 
scribed by the effective Hamiltonian ([1]) using the DMRG 
method^^^^^. Namely, we suggest that the Mott core can 
be employed as a model system to widely study the mag- 
netism in interacting spin models. One of the advantages 
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FIG. 1: (a) The trap-potential strength V dependence of the 
particle density profile ntot{i){= ^"ii^)) ^ completely polar- 
ized fermionic gas with slice pictures for two cases (b) V/t = 4 
and (c) 10. 



using the equivalence is that one can easily reach a very 
high field range. Moreover, although this paper concen- 
trates on the one-dimensional and two-component Fermi 
atom system as a trial problem, higher dimensional, frus- 
trated, and large S cases are also possible to study. 

The starting model Hamiltonian*^'^^ describing trapped 
two-component Fermi atoms under the ID strong opti- 
cal lattice is given by the ID Hubbard model with the 
harmonic trap. 



i^Hubbard = X] {4a^j<T + H.C.^ + ^ 

where, the summation for the pseudo-spin a is taken over 
two-components assigned as a =t and |, respectively. 
cl^ is the creation operator of a Fermi atom with the 
pseudo-spin a at the i-th lattice state, and nicr{= (^la^ia) 
is the site density one for the same pseudo-spin. In the 
first term of the Hamiltonian ([2]), t describes the nearest- 
neighbor hopping parameter and the summation (i, j) is 
taken over the nearest-neighbor sites, and U (> 0) in 
the second term is the on-site repulsive interaction. The 
last term in Eq. ([2|) describes a harmonic trap potential, 
where V is the potential height at the edge sites. N is the 
total number of lattice sites, and Ny is that of fermions 
with cr =t and [ {N^ = + ^^i)- Throughout this 
paper, an atom component with a =t is always a major 
one. As a main numerical method, we employ the DMRG 
to explore the ground state of the model ([2]). At first, 
the number of states kept (m) in DMRG is selected by 
a comparison of the ground state energy with the exact 
diagonalization method for small size {N = 20). In larger 
sizes, we select m which gives no significant difference by 
increasing m further. For N = 60(120) and 180(240) 
in the Hubbard model, we confirm that m = 100 and 
m = 300 is enough, respectively. In addition, for TV = 60 
in the Heisenberg model, m = 100 is selected due to the 
same reason. 

Let us show DMRG results of the model ([2]). Firstly, 
we show atomic density profiles in the case of the per- 
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FIG. 2: The trap-potential strength V dependences of the 
particle density profiles ntot(0(= ''^t(^) + ^i(O) — 66 

and Ni = 4 in (a) U/t = 4 with two slice pictures at (b) 
V/t = 12 and (c) 20 and (d) U/t = 8 with the same ones at 
(e) V/t = 12 and (f) 20. 



feet polarization (P = {N^ - Ni)/N = 1) in Fig. [H 
When V/t > 5, we find the insulating core in the cen- 
ter of the trap, over which the unit filling is spread. 
Since the compressibility is zero and the polarization is 
perfect in this insulating core, it is regarded as a ferro- 
magnetic insulator. We note that this insulating state 
originates from only the Pauli's exclusion principle and 
differs from the Mott state caused by a repulsive inter- 
action between fermions. We also point out that this 
ferromagnetic insulating core can be described by the 
antiferromagnetic Heisenberg model ^ in the presence 
of an infinitely strong magnetic field. 

Next, let us study cases in which the minority spin 
component slightly increases from the zero (the complete 
polarized one). We examine density profiles in two typ- 
ical situations, i.e., those in the presence of relatively 
weak and strong repulsive interaction. The upper and 
the lower panels in Fig. [2] are V/t dependences of density 
profiles of the former (U/t = 4) and the latter (U/t = 8) 
cases, respectively. As seen in Fig. [2l the unity core is 
broken in the central region about above F/t = 10 in the 
weak interaction case (Fig. [2ja)), while its flat plateau 
feature is still kept up to V/t = 20 (Fig. El^d)) in the 
strong interaction one. Here, we note that the break- 
down of the unity core is also observed above V/t = 20 
in the strong interaction {U/t = 8) case. Namely, the 
V-dependent changes in the density profiles are qualita- 
tively equivalent in both the cases. On the other hand, we 
find from these results that the unity core is the so-called 
Mott state since its phase stability actually depends on 
the interaction strength. 
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FIG. 3: The profile changes in the Hubbard model (|2} with 
decreasing the population imbalance ratio at U/t — 8 and 
y/t = 20. In the fixed total particle number A^f = 70, (a) 
= 60, (c) 46, and (e) 38. The spin densities in the 60-sites 
Heisenberg chain with the open boundary condition in the 
external magnetic field given by Eq. ([T]) are plotted with (b) 
gi^BHe^t/J = 0.90, (d) 0.47, (f) 0.15. 



Now, let us concentrate on pseudo-spin structures in- 
side the Mott-phase as seen in Figs. [2fe) and[2](f). Be- 
fore the Mott phase destruction occurs, we find that the 
minority makes a profile like Wigner lattice inside the 
Mott core. The number of the peak in the minority pro- 
file is the same as that of the minority atoms as seen 
in Figs. [2je) and EJf), where the number is just four 
(see Fig. [3l^a) for another case in which the number is 
ten). These Wigner lattice like profiles can be explained 
by the antiferromagnetic Heisenberg model ([T]) in finite 
but strong magnetic field. The effective model ([T]) then 
predicts the spin density wave (SDW) state whose peri- 
odicity is characterized by 2kF = 7r(l — m) where fh is 
the magnetization normalized by the saturated magne- 
tization and kp is the Fermi wave vector in the equiva- 
lent spinless fermion systen>i^. In the present imbalance 
system, since the m is a controlable parameter via the 
population imbalance, the periodicity is given by 



2kf 



where A^,^°** and 



^^Mott _ jyMott' 
^Mott 



(3) 



^Mott 



are the numbers of the 



and down-spin particles participating the Mott core, re- 
spectively, and N^^^^ is the number of the lattice sites 
occupied by the Mott core. Thus, one finds why the 
minority profile shows Wigner crystal like ones, e.g., 
2kF = 7t[1 - (56 - 4)/60] = 27r • 4/60 in Fig. ^f) and 
2kF = 7r[l - (50 - 10)/60] = 27r • 10/60 in Fig. [3l^a), 
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FIG. 4: The comparison of k vs. the Fourier transformed 
spin density ns{k) for Figs. [3ja) andfSlJb), i.e., the trapped 
Hubbard model Eq. ([2} and the Heisenberg model Eq. ([T]) 
with the open boundary condition, in which k = (ni) /(L + 1) 
(i — 1,2,...), where L = 60. In the case of the trapped 
Hubbard model, only the central core region is used for the 
Fourier transformation. 



where the Mott phase covers 60 sites (A^^ott ^ gQ^^ 
^Mott _ _ (^]^Q majority particles contribute to 
make the metallic wings), and TV^^^** = as seen in 
Fig. Wis) Fig. [3l^a). These profiles are really con- 
firmed by the DMRG calculation of 60-sites Heisenberg 
chain model in a magnetic field with the open boundary 
condition e.g., compare Fig. [3]^a) with Fig. [3](b). This 
result clearly demonstrates that the imbalanced Mott 
phases in the trapped Fermi lattice systems are equiva- 
lent with the effective interacting spin model under mag- 
netic field. 

Let us compare the spin density distributions of the 
Mott core with ones of the Heisenberg model in more 
details. For the purpose, we evaluate the Fourier com- 
ponent n^{k) of the spin density distributions n^{i){= 
n^{i) — ^iii)) in a central range (from z = 31 to 90, 
i.e., L = 60) shown in Figs. [3l^a)-[3Kb). Figures [H show 
k = {7r£)/{L + 1) (£ = 1, 2, • • • , L) vs. n^ik). In these 
figures, one can find that a main peak characterizing the 
SDW structure (e.g., i = 21) and other profiles almost 
coincide between both cases. This result indicates that 
the Mott phases confined inside the harmonic trap can 
be well described by the effective Heisenberg model with 
the open boundary condition. 

We further decrease the population imbalance P, i.e., 
increase the number of the minority atoms. Then, the 
results, e.g.. Fig. [3](c) reveals that the SDW periodic- 
ity is reduced according to 2kF = 7r(l — m). We also 
note that by further imbalance decrease, in addition to 
the SDW spin configuration, another modulation struc- 
ture with a wave length being much longer than the lat- 
tice constant appears [see Fig. [3fe)]. This is regarded 
to emerge as a boundary effect since the incommensura- 
tion of 2kF = 7r(l — fh) with the lattice becomes visible, 
i.e., a beating modulation whose periodicity given by 7rm 
is exposed. For example, 2A:f = 7r[l — (33 — 27) /60] = 
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FIG. 5: The trap potential strength V dependences of profiles 
of the particle density ntot{i){= ^^{'^) — 36 

and Ni = 34 in (a) U/t = 4 with two slice pictures at (b) 
V/t = 12 and (c) 20 and (d) U/t = 8 with the same ones at 
(e) V/t = 12 and (f) 20. 



7r(l — 6/60) in Fig. [3l[e), where it is noted that both 
the majority (5 particles) and the minority (5 particles) 
equally contribute to the metallic wing. See another case, 
2kF = 7r[l - (31 - 29)/60] = 7r(l - 2/60) in Fig. [5][f), 
where 5 majority and 5 minority particles also partici- 
pate the metallic wing similar to Fig. ^e). As shown in 
Figs.[3](d) and[3](f), the change of the spin structure seen 
in Figs. [3l^c) and in the Hubbard model ([2]) can be 
also well reproduced by decreasing the strength of the 
magnetic field in the effective model ([1]). One finds that 
even the beating modulation due to the boundary effect 
is also reproduced. 

Let us turn to further small imbalance cases close to 
the balanced one. The upper and lower panels of Fig. [5] 
show V/t dependent profiles in which the population ra- 
tio is 36 : 34 in strong and weak U/t^ respectively. The 
profile in the weak interaction shows that the Mott insu- 
lator core is broken about above V/t = 12 and the almost 
antiferromagnetic staggered profile is lost in the broken 
region as shown in Fig. EJc). The loss of the staggered 
structure is also observed in the periphery^^ around the 
Mott core as seen in Fig. [5fb) (see Fig. [6](a) for another 
case). These results clearly reflect that the staggered 
profile, i.e., the SDW phase is formed only by the spin 
degree of freedom. The staggered profile diminishes in 
the metallic region in which the charge degree of free- 
dom is alive. In addition, inside the Mott core, another 
long modulation is also observed in both the weak and 
strong interaction cases as seen in Figs.[5fb) and[5l[f). In 
order to check the size dependence of this modulation, 
we examine the profiles by simply increasing both the 
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FIG. 6: The site number N and the total atom number A^f 
dependences of the atom profiles with keeping the imbalance 
ratio for (a) = 22 and N = 60, (b) = 44 and N = 120, 
(c) = 66 and N = 180, and (d) 
these cases, U/t — 20, and V/t — Q. 
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lattice sites and the number of total atoms with keeping 
the population imbalance ratio a constant. The mod- 
ulation and its wave periodicity is found to be almost 
size-independent within the range as seen in Figs. El^a)- 
md). These results indicate that such a modulation is 
clearly observable in 1-D atomic Fermi gases loaded on 
optical lattices. In addition, we note that the effective 
Heisenberg model with the open boundary condition can 
reproduce these results. 

We investigated the repulsively-interacting polarized 
1-D Hubbard model with harmonic confinement poten- 
tials by using the DMRG method. Inside the core phase 
(where the site density equals to the unit) emerged uni- 
versally for arbitrary P, we found that its spin structure 
is described by the antiferromagnetic Heisenberg model 
in magnetic field. This equivalence was confirmed by 
DMRG calculations for both the original and effective 
models. We suggest that the repulsively-interacting po- 
larized trapped lattice fermion systems offer various play- 
grounds of not only the Hubbard type but also the in- 
teracting localized-spin one. This idea may have a new 
impact on studies of the magnetism in the solid state 
physics. 
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